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Several analytical mechanical slosh models for a cylindrical tank with flat bottom are 
reviewed. Even though spacecrafts use cylinder shaped tanks, most of those tanks usually 
have elliptical domes. To extend the application of the analytical models for a cylindrical 
tank with elliptical domes, the modified slosh parameter models are proposed in this report 
by mapping an elliptical dome cylindrical tank to a flat top/bottom cylindrical tank while 
maintaining the equivalent liquid volume. For the low Bond number case, the low-g slosh 
models were also studied. Those low-g models can be used for Bond number > 10. The 
current low-g slosh models were also modified to extend their applications for the case that 
liquid height is smaller than the tank radius. All modified slosh models are implemented in 
MATLAB m-functions and are collected in the developed MST (Mechanical Slosh Toolbox). 


Nomenclature 

R q : tank radius 

1i : liquid height 

til n : the n slosh mass 

W n : the n slosh mode frequency 

C, n : the n slosh mode damping 

h n : the n slosh mass location from liquid level 

l n : equivalent pendulum length for n h slosh mass 

Z : the n slosh mass displacement 

f : the coupling force between rigid and slosh dynamics 

f j. : the coupling force between flex and slosh dynamics 

V : characteristic velocity 

p : liquid density 

C T : liquid surface tension 

We : Weber number 

Bo : Bond number 

Ga : Galileo number 
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a : vehicle axial acceleration 

(j ) : velocity potential 

r,0, Z- axes of tank-fixed coordinate system in Figure 3 
X : m zeros of the derivative of the Bessel Function of the 1 ” kind, J { 

t : time 

p : fluid pressure 

V : fluid velocity 

m ljq : liquid mass 

k v : kinematic viscosity of the liquid. 

d : dome depth 

l : total tank length 

R 0 : equivalent tank radius with dome correction 

T * . ... 

n : equivalent liquid height with dome correction 

TJ : wave height above meniscus as shown in Figure 3 

8 : nondimensional wave height above meniscus ( £ = Tj / R 0 ) 

/ 3 : nondimensional meniscus height at tank wall 

a n : expansion coefficient in series for (j) 


b n- 

CK 

C2„ 

C3„ 


expansion coefficient in series for s 


Fourier-Bessel coefficients 


: dimensionless natural frequency of n th slosh mass 


I. Introduction 

When modeling the dynamics for a rocket-propelled spacecraft, it is very important to properly predict the 
behavior of liquid in a propellant tank. There are two approaches to model the liquid sloshing dynamics. One uses 
CFD (Computational Fluid Dynamics) theory, while the other approximates the liquid behavior using a spring-mass- 
damper mechanical model. The former approach is usually used to predict slosh mass motion in time-domain 
simulation. The latter is used in the frequency domain analysis to predict the interaction between slosh and other 
space vehicle dynamics. In this research, several mechanical slosh models ([1] to [4]) are reviewed. The high-g slosh 
parameter model [1] has been derived for a cylindrical tank with a flat top and bottom. Even though spacecraft use 
cylinder shaped tanks, most of those tanks usually have elliptical domes. To extend the application of the analytical 
model in [1] for a cylindrical tank with elliptical domes, modified high-g slosh parameter models are proposed in 
this report by mapping a elliptical dome cylindrical tank to a flat top/bottom cylindrical tank while maintaining the 
equivalent liquid volume. Several analytical low-g slosh parameter models ([2] to [4]) are also studied in this 
research. Under the assumption that the liquid height is larger than the radius of the tank, those low-g models ([2] to 
[4]) were derived for a flat bottom shaped cylindrical tank. In this research, those low-g parameter models are 
modified to extend their application for the case that the liquid height in the tank is smaller than the radius of the 
tank. The same technique which extends the application of the high-g slosh parameters model to account for the 
dome area can be applied to the modified low-g models. The contents of the report are outlined as follows. First, the 
high-g slosh parameters model is briefly reviewed. The modified high-g slosh parameters model is derived next. The 
low-g models for a cylindrical tank with a flat top and bottom will also be reviewed. The necessary modification 
which extends the application of the reviewed low-g models is also introduced. One modification is for the case that 
liquid height is less than the radius of the tanks while the other one is for the tanks with elliptical domes. All 
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modified high-g and low-g slosh models are implemented in MATLAB m-functions and are collected in the 
developed MST (Mechanical Slosh Toolbox). A brief summary of the code flows in terms of the equations to get to 
the slosh parameters are also provided in this report. The outputs of MST are the slosh parameters, which include 
slosh frequency, slosh mass, slosh location and smooth wall slosh damping, as functions of the liquid level in the 
corresponding mechanical slosh model. 


II. Mechanical Slosh Models 

The main interest of the slosh dynamics is the lateral sloshing motion or oscillation of the liquid mass. The 
mechanical slosh model approximates the motion of the movable liquid mass with a rigid mass plus a mass -spring- 
damper model [1] (Figure 1). 


A 


h 


y 

Figure 1. Mechanical Models of Sloshing for a Cylindrical Tank 



2 R, 


The corresponding Laplace transformed mechanical slosh model can be mathematically formulated as [5] 

m n O 2 + 2£ n w n s + w 2 n )Z„ = f rs +f fs ( 1 ) 

The slosh mass, frequency, damping and location in the above equation are functions of physical 
properties of the liquid, liquid level, the tank radius, and the axial body force acting upon the space 
vehicle. 


III. Hydrodynamic Regimes 

The motion of liquid mass was driven by the resultants of the capillary force, body force, inertia force and 
viscous force. The relativity of each force to the others determines the hydrodynamic region [1]. As shown in Figure 
2, the relative importance of inertia and capillary force is defined as Weber number, i.e. 


3 

American Institute of Aeronautics and Astronautics 




Figure 2. Hydrodynamic Regimes [1] 


TT7 Inertia pVR,, 

We = = ( 2 ) 

Capillary a 

Similarly, the acceleration to capillary forces ratio is defined as Bond number, i.e., 

r , Accelerati on paRz 

Bo = = - — - (3) 

Capillary a 

The dimensionless Bond number will determine the shape of the equilibrium surface of the liquid tank as shown in 
Figure 3 [3], For Bond number > 1000, the equilibrium liquid surface is reasonably flat. This scenario is classified as 
a high-g slosh problem. The contact angle is defined as the angle between the water line and the wall. In a high-g 
problem, the contact angle is equal to 90 degrees. As the Bond number decreases, the equilibrium interface starts 
curving. When the Bond number approaches 0, the water line will be parallel to the wall, thus the corresponding 
contact angle is equal to 0 degrees. 

Z 
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IV. High-G Slosh Parameter [1] 


For an incompressible flow, the velocity potential <j)(r,0,z ) must satisfy Laplace’s equation everywhere in 
the liquid volume, i.e. 

vV = o 

and two boundary conditions 

d(j) 


dr 

d(f> 

dz 


= 0 


(4) 

(5) 


)•=«„ 


= 0 


(6) 


Z=~h 


The steady state solution to the above Laplace equation satisfying those boundary conditions is 

cosh(/t m - -^) 

tVn 


</>(r, 0, Z) = X A mn Jn (K COS (" 6 ) , 

0 cosh(/l — ) 

R n 


(7) 


where the values of coefficients A can be solved for a given set of boundary conditions. For an irrotational flow, 
the slosh liquid motion can be described by the unsteady form of Bernoulli’s equation 


+ az(r, &,t) + — + — = 0 

dt 2 /? 


( 8 ) 


The corresponding time-dependent solution is 

O (r,0,z,t) = (/)(r,0,z)e lwt (9) 

Linearizing the unsteady form of Bernoulli’s equation at the free surface and applying the following “Kinematic” 
condition 

d5 50 ll 

= w= for Z — — 

dt dz 2 

the time derivative of the linearized Bernoulli’s equation becomes 

5 2 0 50 „ 

— — + a = U tor 

dt 2 dz 

Plug Eq.(9) into the above equation and solve for w 


h 

z = — 
2 


( 10 ) 


( 11 ) 


= K: TT tanh(A m 


R, 


R n 


( 12 ) 


l o iV o 

The corresponding slosh location ( h ) measured from liquid level and slosh mass ( in ) can be derived as [1] 


2 R () tanh| 


m n = m n q 


'a a ' 

' R„ / 




K =^ptanh| 


A, 


h 


2 R, 


(13) 


(14) 


o J 


The slosh damping due to the physical baffles installed inside the tank can be estimated by using the energy 
dissipation rate of a force free oscillation. For oscillation, the slosh damping can be approximated with [6] 
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<r = 


V 4 n 1 + A 2 


where A is defined as the logarithmic decrement 

A = In 


^ Current Peak Resp onse ' 


(15) 


(16) 


Next Peak Response 

Eq. (15) is usually used to interpret experimental measurement. The analytical slosh damping formulas for anti-slosh 
devices are available in [1] based on fluid dynamics analysis. Liquid slosh in a cylindrical tank without an anti-slosh 
device is damped by the viscous stresses. This type of damping is labeled as the smooth wall damping. The model of 
the smooth wall damping of the 1 st slosh mass for a circular cylindrical tank has been derived [1] 

f h ^ 


h <R 0 Ci = 0 . 19 yfRe 1 + - 


0.318 


sinh 


1.84 


h 


R 


Oj 


1 + 


1 — 

R n 


cosh 


1.84 


h 


R, 


0 J ) 


(17) 


h >R 0 Ci = 0 J 9 jHe 

where the dimensionless parameter Re is defined as 


Re = 




(18) 


The mechanical slosh parameter models cited in the section are valid if slosh displacements are small. If the wave 
amplitude due to the axial load is greater than 10% of the tank diameter, splashing, breaking waves and even rotary 
sloshing will occur, thus the apparent slosh damping will increase in the mechanical slosh model [1]. 


V. High-G Slosh Parameters for a Cylindrical Tank with Elliptical Domes 

The slosh models cited in the previous section are derived for a cylindrical tank with a flat bottom as shown in 
Figure 1. Most propellant tanks in a space vehicle have both elliptical top and bottom domes as shown in Figure 4. 
The slosh parameter models for a flat bottom cylindrical tank cannot be directly applied to the dome area. One way 
to resolve this issue is to map a cylindrical tank with top and bottom dome to a cylindrical tank with flat top and 
bottom while maintaining the equivalent liquid volume [7], Given a tank radius of R 0 , a dome depth of d and a total 
tank length of /, the formula for the dome area correction follows. 
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n 

For the case of liquid level in the top dome area, the equivalent tank radius ( n (| ) and liquid height ( n ) for a 
cylindrical tank with flat top and bottom can be derived as 


R 0 -R 0 . 


(l-h)(2d + h - /) 
d 2 


h* = 


2 d 3 -3 Id 2 + (3 h 2 - 6 hi + 3 l 2 )d + h 3 - 3 h 2 l + 3 hi 2 - l 3 


(19) 


(6 h - 61) d + 3 h 2 - 6 hi + 3/ 2 

r> * 

When the liquid level is in the barrel section as shown in Figure 5, the equivalent tank radius ( R 0 ) and liquid height 
( /? ) for a cylindrical tank with flat top and bottoms can be derived as 


h* = h-- 
3 


( 20 ) 


_ * . * 

For the case of the liquid level in the bottom dome area, the equivalent tank radius ( R 0 ) and liquid height in) for 
a cylindrical tank with a flat bottom can be derived as 


K = R q 


h(2d - h) 


h* = 


d 1 

h(3d -h) 
6 d - 3 h 


(21) 



VI. Low-g Slosh Parameters for a Cylindrical Tank 


Equation (7) in the high-g parameter section can be made nondimensional by replacing r / R 0 with R and 
setting « = 1 ([2]-[4]). 


<p(R,d,Z,r) = ^a m (r)/ 1 (/ i m R)cos(0)e A "' z (22) 

m = 1 


The nondimensional wave shape ( S = T] / R 0 ) can be modeled as [3] 


s(R, 9 , r) = Yjb m (j)J i (A„ R ) cos(0) (23) 

m = 1 
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where T] is the wave height above meniscus as shown in Figure 3. Based on the boundary conditions, the motion of 
the free surface and the fluid velocity at the free surface must be equal to each other, i.e. 

ds dtp dF ( dtp ^ 


dr dZ dR 




-0, 


Z - F 


r 


r £ — 

dr N 


BO 


J _d_ 

R dR 


V 


R 


ds 

dR 


1 


f dF > 


2 A 


V dR j 


K- 


J_A 

R 2 dd 


ds 

Iw 


f dF > 


2 A 


V dR j 


-X 0 Q 2 Rcos OsinQ? = 0, 

F - P 


Z - F 


where 


fi is the largest positive root of the following polynomial 

P 3 N bo -j3 2 - 2/3 = 0 

Solve the two boundary conditions (Equations (24) and (25)) 



( r ^ 

1/2 

1- 

1 T 



l R l) 





n = 1 


m = 1 


X \b n+ tC2 


n = 1 


n = 1 


a + V C3 b 

mn m mn n 

m = 1 


2X Q 2 

/> 1 .. , sin i (K R )cos 0 = O 

1^-1 PMn) J 


and combine the above two equations by eliminating b 


where 


X C2 <«« a m - X C1 - - X C3 ™ X C1 


m=l 


m=l 


j=l 


m = 1 


2Z 0 Q 




cosQf 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 


Clmn — 

C2 , ml = 
C3,„„ - 




x7 1 (2„/?)e^[ 1_(l_s3)1 ' 2 W 

2^ f 1 R 

W-tfp wf Jo A^ BO (l-/? 3 + (9/ 4)/? 2 /? 4 f 2 
X {(l - tf 3 } /2 4 7, a m /?) + (9 / 4)/7 2 /? 2 (1 - tf 3 )7, (47?) - 


9J3 2 R 3 (1 — 0.25/f J )(1 — /? ) 
(l-T? 3 +(9/4)/7 2 /? 4 ) I/2 


3 \l/2 


J[{Z m R)JM„R)}d R 


(31) 


Solve Equation (30) for a n 
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<2 = 


In 


2 n 


s 2 -Cl] s 2 -Cl; s 2 -Q 2 m 

then plug the above solution into Equation (28) and solve it for b ] 

Q\n i Qln - i Q mn 


X 0 Q 3 


b„ = 


2 2 2 /-v 2 

s — £ 2 X s -£2 2 


2 ^2 


Z 0 Q- 


(£> 2 <£> 2 <--<£ 2 ) 


(Qf<Q 2 < -<Q 2 J 


(32) 


(33) 


The slosh parameters in a low-g environment are 

w = Q I — 

" n \l R 


0 


OT„ = — 


M 


/„ - 




N 


H 


O 2 

bo y 


(^o) 


M 


(34) 


Y-^JXX )e KP , H =YO JX / 1 ) 

* 1 V m / ’ ft *~-nm IV fti 7 

m=l 


W=1 


»,(A- 


1 


n 




1V1 

I 

/ft=l 




m / 


P nJ\(K ) e 

X. 


x m p 


-fi 


P~ 


N 


H„ 


o 2 

BO J 


/?n 


Notice that the above equation predicts the slosh motions in a cylindrical tank when the liquid height is larger than 
the tank radius (i.e. h > R„). In the next section, modified slosh model is proposed to extend its application for the 
case of h < R„ 

The high-g smooth wall damping model (Eq. (17)) cannot used to predict slosh damping in a low-g environment. 
The low-g smooth wall damping of the 1 st slosh mass for a circular cylindrical tank has been provided [4] 

[0.83(jV o J-'' 2 +0.096(Ar JO )~' n , N m >10 
\ 0.83 (/V m ) "(l + 8.21 :V„„ ) "4 otherwise 
where the dimensionless parameter Galileo number Nqa is defined as 


ft 


(35) 


n ga = 


Qi/?o 


(36) 


(0.4647 ) 2 k v 

VII. Low-G Slosh Parameters a Cylindrical Tank with h/R 0 Correction 

For the case that the height of liquid level is less than the radius of the tank, the slosh parameter models 
derived by Dodge [3] must be modified as follows [7], 



OT„ = — 




(npR\)tanh(^—) 


R n 


(37) 


1 npRl 


R 0 n; 


m„ 


M 

I 

m=l 


1 


A 


P~ T 

A m J 


P nJMm ) e 

2 




-p 


p- 


N 


H 


BO J 


Q: 


A h 
>R 0 tanh ( — ) 
2 R n 


Notice that Equation (37) predicts the slosh motions in the cylindrical tank with a flat bottom. For a cylindrical tank 
with elliptical top and bottom domes as shown in Figure 4, the tank radius and liquid height in the slosh parameters 
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must be modified using Equations (19), (20) and (21) for the liquid height in the upper dome, barrel, and lower 
dome, respectively [7], 

The derivation of slosh parameters in Equation (34) assumes that the contact angle is zero based on the linearized 
approximation [3], For the 90 degree contact angle case [1], the low-g slosh frequency is 


2 

W„ = 


2 nl & 

n A 7 

a 

+ — 

A n tanh 

r 

p R l 



, r 0 . 


(38) 


If the acceleration a reverses the direction such that B a < — ) , then the frequency is imaginary which means 


the interface will be become unstable. The corresponding critical acceleration 

f _ '\ 


a. 


=w : 


(39) 


P R o J 

is the minimum adverse acceleration to maintain surface stability under zero-g environment [1], 

Unlike Eq. (17), a correction equation involving h/R 0 has not been attempted for low-g slosh damping model (Eq. 
(35)) due to insufficient data collected [4], 


VIII. Mechanical Slosh Toolbox (MST) 

All mechanical slosh models derived in the previous sections have been implemented in MATLAB m- 
functions and are collected in the MST (Mechanical Slosh Toolbox). The assumptions for each model are 
summarized in the Table 1. The MST can be used to predict slosh dynamics in various body forces (low-g or high-g) 
by assuming a propellant tank is cylindrical with upper and lower domes. The output of MST are the slosh 
parameters, which consists of slosh frequency, slosh mass, slosh location and smooth wall slosh damping as 
functions of liquid fill level, in the mechanical slosh model. Two scripts are provided to help end users to generate 
slosh parameter data. 

• mst_slosh_table.m : generates low-g/high-g slosh tables 

• mst_slosh_driver.m : plots low-g/high-g slosh parameters 


Model 

High G in [1] 

Low G in [3] 

High G in MST [7] 

Low G in MST [7] 

Bond Number (Bo) 

Bo > 1000 

10 < Bo < co 

Bo > 1000 

10 < Bo < co 

Cylindrical Tank End 
Dome Assumption 

No Domes 

No Domes 

Rounded Domes 

Rounded Domes 

Liquid Height Limit 

None 

h>R 0 

None 

None 

Contact Angle 

90 degrees 

0 degrees 

90 degrees 

0 degrees 

Slosh Parameters 

Mass, Frequency, Mass 
Location, Smooth Wall 
Damping 

Mass, Frequency, Mass 
Location, Smooth Wall 
Damping 

Mass, Frequency, Mass 
Location, Smooth Wall 
Damping 

Mass, Frequency, Mass 
Location, Smooth Wall 
Damping 


Table 1. Assumptions on Mechanical Slosh Models 


To demonstrate the uses of the MST, a cylindrical tank with liquid Oxygen propellants filled is assumed. The tank 
has both elliptical top and bottom domes with arbitrary tank radius of 8.8 feet and dome height of 9.8 feet. For a 
high-g flight environment, the “mst slosh table.m” generates a 1-D lookup table as shown in Figure 6. The lookup 
variable is LiqLevel. In this example, 20 difference liquid levels are specified between 0.1 ft above the bottom of the 
tank and 0. 1 ft below the top of the tank. 
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LiqLevel (ft) 

Freq(Hz) 

Mass (slug) 

Zeta (-) 

hi (ft) 

1 . 00000e-001 

2 . 957 65e-001 

4 . 59062e-001 

3 . 58209e-003 

5. 00634e-002 

2 . 74285e+000 

3 . 17044e-001 

2 . 96167e+002 

2 . 53919e-004 

1 . 42377e+000 

5 . 38571e+000 

3 . 41341e-001 

9 . 39358e+002 

1 .54037e-004 

2. 91349e+000 

8 . 02856e+000 

3. 68245e-001 

1 . 60573e+003 

1 .17351e-004 

4 . 53824e+000 

1 . 06714e+001 

3. 94811e-001 

1 . 96746e+003 

1 .00994e-004 

6 . 21771e+000 

1 .33143e+001 

4. 06851e-001 

2 . 08929e+003 

9 . 44904e-005 

7 . 48280e+000 

1 . 59571e+001 

4. 10912e-001 

2 . 13121e+003 

9 . 15112e-005 

8 . 30691e+000 

1 . 86000e+001 

4. 12265e-001 

2 . 14526e+003 

8 . 99413e-005 

8. 81813e+000 

2 . 12428e+001 

4. 12713e-001 

2 . 14993e+003 

8 . 78210e-005 

9 . 12568e+000 

2 . 38857e+001 

4. 12862e-001 

2 . 15148e+003 

8 . 78210e-005 

9 . 30730e+000 

2 . 65285e+001 

4. 12911e-001 

2 . 15199e+003 

8 . 78210e-005 

9 . 41339e+000 

2 . 91714e+001 

4. 12927e-001 

2 . 15216e+003 

8 . 78210e-005 

9. 47496e+000 

3 . 18142e+001 

4. 12933e-001 

2 . 15222e+003 

8 . 78210e-005 

9 . 51056e+000 

3 . 44571e+001 

4. 12935e-001 

2 . 15224e+003 

8 . 78210e-005 

9 . 53110e+000 

3 . 70999e+001 

4. 12935e-001 

2 . 15224e+003 

8 . 78210e-005 

9 . 54293e+000 

3 . 97428e+001 

4. 12935e-001 

2 . 15225e+003 

8 . 78210e-005 

9 . 54975e+000 

4 . 23856e+001 

4. 16244e-001 

2 . 05161e+003 

8 . 88786e-005 

9. 40410e+000 

4 . 50285e+001 

4 . 36733e-001 

1 .53776e+003 

9 . 55211e-005 

8 . 54553e+000 

4 . 76713e+001 

4. 95292e-001 

7 . 22800e+002 

1 .15363e-004 

6. 64437e+000 

5 . 03142e+001 

1. 09276e+000 

6 . 26656e+000 

3 . 78063e-004 

1 . 36497e+000 


Figure 6. High-g Table from “mstsloshtable.m” 


For a low-g flow flight scenario, the “mst slosh table.m” generates a 2-D lookup table as shown in Figure 7. The 
lookup variables are LiqLevel and Bond Number (or Gravity). In this example, 6 difference liquid levels are 
specified between 0.1 ft above the bottom of the tank and 0.1 ft below the top of the tank. In addition to the liquid 
level variable, 6 logarithmically equally spaced Bond numbers are also assigned between decades 10 1 and 10 6 to 
complete the 2-D lookup variables. 
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LiqLevel (ft) 1 . OOOOOe-OOl 1.01428e+001 2.01857e+001 3.02285e+001 4.02714e+001 5.03142e+001 

Bond Number 1.00000e+001 1 . 00000e+002 1.00000e+003 1.00000e+004 1.00000e+005 1.00000e+006 

Gravity (ft/s A 2) 5.28047e-005 5.28047e-004 5.28047e-003 5.28047e-002 5.28047e-001 5.28047e+000 

Freq (Hz) 

3 . 8 961 6e-004 1.19336e-003 3.76583e-003 1.19414e-002 3.78288e-002 1.19729e-001 

5 . 14334e-004 1.57536e-003 4 . 97129e-003 1.57639e-002 4 . 99379e-002 1.58055e-001 

5 . 43511e-004 1 . 66473e-003 5.25331e-003 1 . 66581e-002 5.27708e-002 1.67021e-001 

5 . 43959e-004 1.66610e-003 5.25764e-003 1.66719e-002 5.28143e-002 1.67159e-001 

5 . 43966e-004 1.66612e-003 5.25770e-003 1 . 66721e-002 5.28150e-002 1 . 67161e-001 

1.43951e-003 4.40911e-003 1.39136e-002 4.41198e-002 1.39766e-001 4.42363e-001 

Mass (slug) 

5 . 12120e-001 4.29155e-001 4 . 41741e-001 4.54312e-001 4 . 60808e-001 4 . 63901e-001 

2 . 14654e+003 1.79879e+003 1 . 85155e+003 1 . 90424e+003 1.93147e+003 1.94443e+003 

2.39699e+003 2.00867e+003 2 . 06758e+003 2.12641e+003 2.15682e+003 2.17130e+003 

2 . 40094e+003 2.01198e+003 2.07099e+003 2.12992e+003 2.16038e+003 2.17488e+003 

2 . 40100e+003 2 . 01203e+003 2.07104e+003 2.12997e+003 2.16043e+003 2.17493e+003 

6 . 99085e+000 5.85830e+000 6.03012e+000 6.20172e+000 6.29040e+000 6.33262e+000 

Zeta ( - ) 

3.34844e-002 1 . 58530e-002 6.55581e-003 2.93673e-003 1.41420e-003 7.20275e-004 

7 . 7781 6e -00 3 1 . 10525e-002 3.85346e-003 1 . 41918e-003 5.61566e-004 2.41014e-004 

7 . 7781 6e -00 3 1 . 10525e-002 3.85346e-003 1 . 41918e-003 5.61566e-004 2.41014e-004 

7 . 7781 6e-003 1 . 10525e-002 3.85346e-003 1 . 41918e-003 5.61566e-004 2.41014e-004 

7 . 7781 6e -00 3 1 . 10525e-002 3.85346e-003 1 . 41918e-003 5.61566e-004 2.41014e-004 

3.34844e-002 1 . 58530e-002 6.55581e-003 2.93673e-003 1.41420e-003 7.20275e-004 

hi (ft) 

4 . 20068e-002 5.05826e-002 5.05170e-002 5.02388e-002 5.01310e-002 5.00913e-002 

4.95198e+000 5.96294e+000 5.95521e+000 5.92242e+000 5.90971e+000 5.90502e+000 

7.56963e+000 9.11499e+000 9.10317e+000 9.05304e+000 9.03361e+000 9.02645e+000 

7 . 96411e+000 9.59001e+000 9.57757e+000 9.52484e+000 9.50440e+000 9.49686e+000 

8 . 0137 4e+0 00 9.64977e+000 9.63725e+000 9.58419e+000 9.56362e+000 9.55603e+000 

1.14531e+000 1.37913e+000 1.37734e+000 1.36976e+000 1.36682e+000 1.36573e+000 

Figure 7. Low-g Table from “mst slosh table.m” 

The “mst slosh driver.m” generates high-g verse low-g comparison plot for the selected tank as shown in Figure 8 
where the high-g slosh parameters are plotted in a blue color and the low-g slosh parameters are plotted in a green 
color. 
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Figure 8. High-g vs Low-g plots “mst slosh driver.m” 


The high level code flows for both scripts are drawn in Figure 9 and Figure 10. The function list for MST follows. 

• any_g_calc.m : Main function used to compute slosh parameters 

• low_g_calc.m : Computes slosh parameters using low-G model 

• low_g_damp.m : Low-G damping model 

• low_g_freq.m : Low-G frequency, mass, and location model (dimensionless) 

• fourier_bessel_coeffs.m : Compute Fourier-Bessel coefficients 

• bessel_first.m : Bessel function of first kind 
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• betaFromNbo.m : Solves for beta 

• high_g_calc.m : Computes slosh parameters using high-G model 

• mst_constants.m : Constants used in the MST 


mst slosh table 



Figure 9. Code Flows for “mstsloshtable.m” 
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mst slosh driver 



Plot Low-g Slosh Paramamters 
Vs. High-g Slosh Paramamters 
(Mass, Location, Frequency, Smooth Wall Damping) 


Figure 10. Code Flows for “mst slosh driver.m” 

IX. Summary 

The mechanical slosh models for a cylindrical tank with a flat bottom have been reviewed under this study. The 
high-g slosh models (Bond number > 1000) from [1] were modified in this study for a cylindrical tank with elliptical 
domes. For the low Bond number case, the low-g slosh models cited in the literature ([2] to [4]) were also studied. 
Those low-g models can be used for Bond number > 10; an additional model study may be needed for the smaller 
Bond number scenario. The low-g slosh models from [2] to [4] were also modified to extend their applications for 
the case that liquid height is smaller than the tank radius and the case of a cylindrical tank with elliptical domes. 
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